Design of UHECR telescope with 1 arcmin resolution 
and 50 degree field of view 
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Abstract 

A new telescope design based on Baker-Nunn optics is proposed for observation of ultra high-energy cosmic rays 
(UHECRs). The optical system has an image resolution of smaller than 0.02° within a wide field of view of 50° 
angular diameter. When combined with a high-quality imaging device, the proposed design enables the directions 
of UHECRs and high-energy neutrinos to be determined with accuracy of better than 1 arcmin. The outstanding 
resolution of this telescope allows charge-separated cosmic-rays to be resolved and the source to be determined 
accurately. This marked improvement in angular resolution will allow the multi-wavelength and "multi-particle" 
observations of astronomical objects through collaboration with established astronomical observations. 

Key words: ultra high-energy cosmic ray; fluorescence detection; the Baker-Nunn optics; 1 arcmin angular 
resolution, wide field of view 



1 Introduction 



The observation of ultra high-energy cosmic rays 
(UHECRs) and high-energy neutrinos (HE-v's) 
with directional accuracy of < 1 arcmin (') is ex- 
pected to lead to revolutionary progress in par- 
ticle astrophysics through "multi-particle" obser- 
vations in cooperation with multi-wavelength as- 
tronomical observations. An excellent example is 
the recent prompt observation of X-rays from a 
gamma-ray burst (GRB) with directional precision 
of 1' by Beppo-SAX, an event that marked the 
start of multi-wavelength observations using many 
advanced astronomical telescopes [1-3]. The di- 
rectional accuracy of 1' is several tens times better 
than that of ~ 1° currently operating, constructing, 
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or planning UHECR detectors [4-7] aim for. 



Upon entering the atmosphere, UHECRs produce 
an air shower (AS) through interaction with air 
nuclei. The AS axis coincides with the UHECR 
incident angle, and the intrinsic angular resolution 
of the shower is 0.3' or smaller for a UHECR of 
10 17 eV. This narrow dispersion results from the 
high Lorenz factor of secondary particles (~ 10 4 ) 
in the direction of the parent particle, and the very 
low transverse momenta of secondary particles 
(rarely exceeding 1 GeV/c) as reported by the AT- 
LAS group [8]. For UHECRs, more precise recon- 
struction of the primary direction can be expected 
due to the higher Lorenz factor. 

ASs can be detected via the isotropic fluorescence 
emission of the air nitrogen that it excites. An op- 
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tical system with imaging devices can be used to 
track the longitudinal development of an AS. This 
technique was first employed in the Fly's Eye de- 
tector [9] and will be part of several future projects 
on UHECR detection. The key advantage of flu- 
orescence detection is the capacity for achieving 
fine angular resolution. The signal to noise ratio 
(S/N) for each pixel of the imaging devices is in- 
versely proportional to the square root of the an- 
gular diameter of the field of view (FOV). A better 
S/N will allow the device to trigger charge lev- 
els, making it possible to analyze less energetic or 
more distant events. 

Galactic and extragalactic magnetic fields are con- 
sidered to have intensities on the order of /iG, nG, 
and below [10]. Thus, UHECRs are expected to 
experience only slight deflection due to galactic 
and extragalactic magnetic fields, and cannot be 
confined by the galactic magnetic field. The ex- 
istence of an "ankle" at E ~ 10 19 eV is often in- 
terpreted as a crossover from a steeper galactic 
component to a harder component of extragalactic 
origin. 

As the possible sources of astrophysical accel- 
erators capable of producing UHECRs are very 
limited, UHECRs that have been deflected only 
slightly can be traced back to their sources. Thus, 
high-resolution measurements of the arrival direc- 
tion will lead to significant progress in cosmic-ray 
astronomy. The magnetic deflection of UHECRs 
during propagation induces charge separation. One 
of the important implications of 1 ' angular reso- 
lution is that this separation of charge can be re- 
solved, allowing unambiguous particle identifica- 
tion to be performed based on the deflection from 
sources. Such observations may also allow us to 
study the structure of extragalactic magnetic fields. 

Based on the above considerations, we propose a 
telescope design that features a FOV of 50 ° x 50 
° and 1 ' resolution for the arrival direction. This 
is realized through a combination of a large-area 
image intensifier (II) connected to a CCD cam- 



era [11] and a customized Baker-Nunn optical sys- 
tem. In this paper, we describe the conceptual de- 
sign of the modified Baker-Nunn optical system. 
The expected performance of the telescopes for 
UHECR detection is also presented assuming that 
the image resolution of the CCD camera is better 
than that of the optical system. 



2 General Concept of Optical System 

Our design requirements can be summarized as 
follows: 

(1) a large FOV of 50° x 50° to cover a large 
solid angle for UHECR detection; 

(2) an image resolution of ~0.02° to achieve 1 ' 
resolution for the arrival direction; 

(3) a large caliber with a diameter of ~2000 mm 
to distinguish faint fluorescence. 

In addition, cost performance will be a critical is- 
sue with respect to construction of a large-aperture 
array of telescopes. 

Among several approaches to meet our require- 
ments, spherical-mirror optics with an aperture di- 
aphragm on the mirror's curvature center appears 
to offer some advantages for the proposed tele- 
scope. This simple optics produces homogeneous 
images within the entire FOV. Typical off-axis 
aberrations, principally coma aberration, are elim- 
inated, although spherical aberration is a factor. 
Thus, from a review of the variants of this optics, 
the Schmidt, Makstov, Baker Super-Schmidt, and 
Baker-Nunn optics were selected as candidates. 
The Schmidt optics has a one-element corrector 
lens at the diaphragm, and provides higher resolu- 
tion than optics without the corrector lens. How- 
ever, we found that it was difficult to achieve the 
required image resolution over the wide FOV. The 
Makstov and Baker Super-Schmidt optical sys- 
tems require thick corrector lenses, and as such 
are not suited to large-scale telescopes. 
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The Baker-Nunn optics has a three-element cor- 
rector lens consisting of convex and concave 
lenses with spherical or aspherical surfaces [12]. 
The basic premise of this optics is that each ele- 
ment has easy-to-make spherical surfaces, residual 
spherical aberration after assembly can be can- 
celed by adjusting the spacing between lens ele- 
ments, and chromatic aberration can be eliminated 
by selecting appropriate lens materials. Figure 1 
shows a schematic of the Baker-Nunn optics. The 
second lens element is an aperture diaphragm. 
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Fig. 1. Schematic of proposed telescope based on 
Baker-Nunn optics. The second lens element (3,4) is 
an aperture diaphragm. 

Although the advantages of the Baker-Nunn op- 
tics are not related directly to the specific area of 
improvement pursued in the proposed design, the 
three-element corrector lens in this optics provides 
a large degree of freedom for optimization. As a 
result, it is considered the most suitable for our 
purposes. In our design, the first and third lens el- 
ements are 2400 mm in diameter, and the second 
element is 2000 mm in diameter, and all elements 
are 10 mm thick. Considering the cost, and com- 
plexity of fabrication processes, we designed the 
corrector lens as an acrylic Fresnel lens as appro- 
priate for a large-aperture telescope. As described 
in the following section, the surface shape of each 
lens element was determined through an optimiza- 
tion process. 



3 Optimization by Ray Tracing 



3.1 Ray Tracing 



We developed ray-tracing program written in C++ 
in order to optimize the optical system and ex- 
amine its performance. The program models the 
refraction and reflection of incident rays on each 
surface, and allows a bundle of on- and off-axis 
rays to be traced numerically. As is common in 
designing an optical system, the shape of the front 
surface of each optical component is expressed as 



h 2 /R 



\ + ^l-h 2 /R 2 
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(1) 



where Z, h, and R represent the surface sag, radial 
height, and curvature of the surface. Z is positive 
in the direction of on-axis incident rays. A, B, and 
C represent the 4th, 6th, and 8th order aspheric 
coefficients. The surfaces of each lens element and 
the mirror are numbered 1 to 7 as shown in Fig. 1. 

The focal surface is defined as the surface mini- 
mizing the image spot for a bundle of parallel rays. 
The spot size is defined as the root mean square 
(RMS) of positions of parallel rays at the focal 
surface. The fineness of the bundle of rays ana- 
lyzed is determined according to the purpose of 
ray tracing: a rather rough interval between rays 
is selected for optimization, and a fine interval is 
selected for performance evaluation. The spot size 
(s) is a function of the incident angle (0) and the 
wavelength (A) of parallel rays. 



3.2 Optimization 



We used the modified Powell's method as the al- 
gorithm for optimization based on the following 
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evaluation function: 

F = £(J>(e,A) 2 ) 2 , 
e x 



where varies from 0° to 25°, and A varies 
from 330 to 410 nm, taking into account the air- 
fluorescence emission spectrum. We adopted the 
fourth power of spot size as the evaluation func- 
tion in order to obtain flat resolution over the 
entire FOV. 

In principle, the curvature, aspheric coefficients, 
thickness, and position of each lens element are 
free parameters in optimization. In order to reduce 
the number of free parameters, we imposed the 
symmetry described in Ref. [12] in which the cor- 
rector lens is symmetric with respect to the plane 
including the center of the second element and 
which is perpendicular to the axis of the ray. In ad- 
dition, the thickness and position of each surface 
were fixed. These constraints allowed much more 
efficient optimization. The optimized parameters 
of the optical system are given in Table 1. The axis 
and origin of the Z coordinate are coincident with 
the axis of the rays and the center of the correc- 
tor lens. Note that we took the effect of the region 
obscured by the focal surface into account in the 
optimization. 

Figure 2 shows the surface sag of lenses 2, and 
3 as a function of radial height. The surface sag 
remained well below 0.5 mm when the width of 
Fresnel grooves was chosen to be 10 mm. Thus, the 
10-mm thickness of each lens element is sufficient, 
although the thickness may be changed according 
to the support method employed for the corrector 
lens. We confirmed that the same level of spot- size 
resolution can be obtained by another optimization 
using a different lens thickness. 

Figure 3 shows the spot diagrams of the opti- 
mized optical system. The shape and size of spot 
diagrams changed with the wavelength and the 
incident angle of the incoming rays. The four 



Table 1 

Optimized parameters of the optical system. 



Surface Z 




Parameters 


1: -405.0 


R: 


oo 


2: -395.0 


R: 


-14767.2 




A: 


2.30448 x 10~ n 




B: 


-7.03756 x 10- 18 


3: -5.0 


R: 


-30479.6 




A: 


1.45344 x 10- 11 




B: 


-6.37472 x 10- 18 




C: 


-8.78219 x 10~ 25 


4: 5.0 


R: 


30479.6 




A: 


-1.45344 x 10- 11 




B: 


6.37472 x 10~ 18 




C: 


8.78219 x 10~ 25 


5: 395.0 


R: 


14767.2 




A: 


-2.30448 x 10- 11 




B: 


7.03756 x 10~ 18 


6: 405.0 


R: 


OO 


7: 3082.0 


R: 


-3060.0 



(unit: mm) 
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Fig. 2. Sag of Surfaces 2 and 3 as a function of radial 
height. The lens surface is constructed from a series of 
concentric grooves (Fresnel lens), each 10-mm thick. 
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Fig. 3. Spot diagrams for optimized optical system. 



strongest fluorescence peaks, A = 336, 358, 391, 
and 425 nm, which dominate the filtered fluo- 
rescence spectrum after propagation through air, 
were selected as wavelengths for the diagrams. 
The spot diagrams for incident angles of 5°, 15°, 
and 25° are shown. 

Figure 4 shows the spot size distribution as a func- 
tion of incident angle, taking the characteristic 
wavelengths in the fluorescence spectrum into ac- 
count. Solid and dashed curves represent the spot 
size optimized to cover a FOV of ±25° and ±20°, 
respectively. The lens parameters shown in Table 1 
are the result of the former optimization, and give 
a spot size of smaller than 0.02° over a FOV of 
±25°. The results were confirmed by cross checks 
using Mathematica Optica [13] with correction for 
the refractive index of the atmosphere. 

In parallel with the development of this optical 
system, the present authors have also developed 
a high-density CCD imaging detector with inte- 
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Fig. 4. RMS distribution of spot size as a function of 
incident angle. 

grated large-area image intensifier, which has been 
demonstrated to have a resolution of ~ 0.02° [1 1]. 
This imaging system is mounted at the focal sur- 
face of the optimized optics to give the proposed 
telescope design. 
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4 Performance 



A solid angle of nearly 2% sr can be observed using 
an array of 8 telescopes (i.e., a station), with 7 tele- 
scopes arranged azimuthally and 1 telescope di- 
rected upward. It should be noted that the telescope 
directed upward is primarily important in HE-v 
detection rather than UHECR detection [14]. 

A Monte Carlo simulation was performed to ex- 
amine the resolution obtained for the arrival direc- 
tion of UHECRs. The event generator calculated 
the signal strength produced by atmospheric fluo- 
rescence emission from ASs with a given geome- 
try, and built simulation events by calculating the 
number of photoelectrons detected by each pixel 
of a CCD camera at each mirror. The simulation 
includes the longitudinal and lateral profiles of the 
AS, the model of the atmosphere, the air fluo- 
rescence yield, and the detector optics. We also 
took into account fluctuation of the first interac- 
tion depth, impact point, and directional angles of 
shower cores assuming appropriate distributions. 

Current and planned fluorescence detectors [5,7] 
have a pixel resolution of ~1°. The angular res- 
olution of the reconstructed track was evaluated 
by analyzing event samples that triggered the cur- 
rent 1° -pixel detectors. Such signals are triggered 
by having at least 6 pixels (1° resolution) fire at 
4a greater than the night sky background, with a 
corresponding reconstructed track extending more 
than 5° in the FOV [5,7]. A night sky background 
level of 100 photoelectrons//! s/(l° x 1°) was as- 
sumed. Figure 5 shows the comparison of event 
displays between the 0.02° and 1° systems. 

Applying a stereo reconstruction of the AS, the ar- 
rival direction is determined by the geometrical in- 
tersection of shower-detector planes (SDPs). Since 
the resolution of the normal vector of a SDP is the 
primary factor affecting the resolution of arrival 
direction, we discuss here for simplicity the accu- 
racy of SDP measurement. The angular resolution 
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Fig. 5. Two examples of AS events for 10 20 eV protons, 
showing events with impact parameter of R p = 53.5 
and 10.9 km. 

of SDP is plotted in Fig. 6 for 10 19 and 10 20 eV 
protons. The angular resolution is determined from 
the integral of the distribution encompassing 68 % 
of events. Figure 7 shows the variation in the an- 
gular resolution of the SDP with proton energy. 
Resolution of < 1' is obtained at above 10 18 eV 
for the SDP. In the case of stereo reconstruction 
with two or more stations, the pointing resolution 
is also expected to be ~ 1' after reconstruction of 
the geometrical intersection of SDPs. 

The S/N of each pixel is proportional to (A/A9)? 
[9], where A represents the light collection area, 
corresponding to the reflection area for simple mir- 
ror optics and the area of the aperture diaphragm 
for the Baker-Nunn optics, and AO denotes the an- 
gular size of each pixel. Note that this proportion- 
ality holds only when the effects of lateral distri- 
bution are neglected. Assuming the same size of A, 
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Fig. 6. The angular resolution of SDP for (Top) 10 
and (Bottom) 10 20 eV protons. 



the maximum values are shown in Fig. 8. Note 
that the event selection applied here differs from 
that applied to compare the angular resolution as 
it is not necessary to apply the 6 pixel condition 
to compare the S/N of the 1°- and 0.02°-pixel de- 
tectors for distant ASs. The solid curve shows the 
expected S/N value for 1° pixels calculated from 
that of 0.02° pixels by the theoretical proportion- 
ality. The expected and obtained S/N values are 
in good agreement for the distant events (R R > 
70 km). For the closer events, the lateral distribu- 
tion of the shower-core image becomes significant 
compared with that for the 1° pixel image (see the 
bottom panel of Fig. 5). As a result, the expected 
S/N ratio is appreciably lower than the simulation 
result. 
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Fig. 7. The energy dependence of the angular resolution 
for SDP determination. 

the S/N value for the 0.02°-pixel imaging device 
is expected to be about 7 times higher than that 
from the 1° -pixel one. However, this is expected 
to be affected by the quantized number of photo- 
electrons at these small dimensions. To examine 
this, the AS profile was simulated by the Monte 
Cairo method and the quantitative improvement in 
the S/N was examined for the 0.02° pixels. The 
AS events of protons with energy of 10 20 eV were 
generated isotropically for a given impact param- 
eter (i? p ). The maximum value of S/N (obtained 
event-by-event) was averaged for a given i? p , and 
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Fig. 8. Averaged maximum values of S/N for 0.02° 
and 1° pixels as a function of 7? p for 10 20 eV protons. 

We also define another signal to noise ratio 
(S/N TRK ) appropriate for the track captured by 
the imaging devices. The signal (5) is the sum of 
charges detected by all pixels associated with the 
track. The background (BG) is the sum of back- 
ground for all pixels associated with the track. 
Then, S/N TRK is defined as S/N TRK = S/\/BG. 
Figure 9 shows the distributions of S/N TRK for 
protons of 10 19 and 10 20 eV. The 0.02° device ex- 
hibited excellent S/N TRK values. The significant 
improvement in S/N TRK was obtained because the 
lateral extent of the AS is smaller than the pixel 
resolution of 1°. Thus, we have demonstrated 
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the improved resolution of the arrival direction 
offered by the proposed system. The detailed per- 
formance of the detector, including aperture, en- 
ergy resolution, and shower-max resolution, will 
be discussed elsewhere. 



T. Aoki for continuous support. We also thank N. 
Manago and M. Jobashi for their useful help in the 
development of the Monte Carlo simulation. 
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5 Conclusion 



We have presented the conceptual design of a new 
type of telescope for UHECR detection based on 
the Baker-Nunn optical system. The telescope has 
a spot size resolution of smaller than 0.02° within 
the entire FOV of 50°, and has been shown to be 
capable of 1 ' angular resolution for the arrival di- 
rection of UHECRs. This design has the potential 
to open new possibilities in the field of UHECR 
and HE-v physics through high-resolution obser- 
vation. 
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